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<^ I present ongoing developments on hybrid models which combine the potential model, 

\q • for the description of the background, with the R-matrix model, for the description of 
discrete states. Physical parameters of discrete states are converted into formal R-matrix 
parameters, taking into account the potential background. The nuclear part of this po- 
tential is constructed by a new variant of an inversion method based on supersymmetric 
quantum mechanics. The method is illustrated on 12 C + a S-wave elastic phase shifts. 
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1. MOTIVATIONS 

In nuclear astrophysics, the key issue is the determination of cross sections at very low 
energies, representative of stellar conditions. Because of the Coulomb repulsion between 
nuclei, these cross sections are sometimes minuscule, in particular for non-resonant pro- 
cesses, and hence very hard to measure. Nuclear theory has thus an important role to play 
for the extrapolation of measured data down to astrophysical energies. A reliable extrap- 
olation requires (i) a precise fit of all available data, (ii) a good account for both resonant 
^ ■ and non-resonant processes, (iii) a small number of parameters to limit the extrapolation 
uncertainty. The R-matrix phenomenological model PQ is the only one to date which is 
able to precisely fit experimental data; however, whereas its description of discrete (res- 
onant) states is satisfactory, its description of the non-resonant background, based on 
phenomenological high-energy wide states, requires a lot of additional parameters and 
lacks physical motivation. 

I demonstrate here the interest of the potential model for a better description of 
the background, complemented by an R-matrix phenomenological description of discrete 
states. Such hybrid models have already been used in previous works |2|3| but the model 
developed here presents two new features. First (section 2), discrete states are character- 
ized by physical parameters (energy and width) rather than by formal R-matrix parame- 
ters, a procedure recently advocated for in the traditional R-matrix method |4I5| . Second 
(section 3), the potential model is constructed by inversion of scattering data, which leads 
to a better- quality fit than phenomenological potentials like the Woods-Saxon one. 

As an example, I apply the method to the S'-wave phase shifts deduced from a 
recent high-precision measurement of the 12 C + a elastic scattering jjj. This analysis 
is preliminary to the analysis of the P and D waves for this system, which dominate 
the capture reaction at astrophysical energies. For these waves, the present method 
could lead to a better description of the interference effects between the background 
and the subthreshold states, which are known to strongly influence the reaction rate 
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at astrophysical energies. The 2 + subthreshold state has recently been studied in the 
framework of the inversion potential model [S] but a generalization to a hybrid model is 
necessary to study the 1~ subthreshold state and to get a precise fit of all available data. 

2. HYBRID MODEL WITH PHYSICAL PARAMETERS 

I use the hybrid potential/R- matrix model presented in refs. |2|3j . The general principle 
is to calculate, for a given partial wave, the R-matrix corresponding to the potential model, 
Rpot(E), and to add to it phenomenological poles corresponding to discrete states. This 
leads to the total R-matrix 



R(E) = R pot (E) + R Tes (E) = R pot (E) + ]T -rr^, (1) 



where E is the center-of-mass energy and only one channel is considered for simplicity 
(this is sufficient for 12 C + a elastic scattering). Discrete states are characterized by their 
formal energy E\ and reduced width 7 2 . The main advantage of the above formula is its 
simplicity, which also implies a very convenient use in numerical calculations. Its drawback 
is that the formal parameters cannot be directly interpreted as physical state properties: 
their value depend on the R-matrix radius and boundary parameter B. This creates strong 
parameter correlations in R-matrix numerical fits where the R-matrix radius itself can be 
varied. It is thus preferable to use physical parameters E\, 7 2 as fitting parameters and to 
convert them into formal parameters for the actual calculation of the fitted quantity. In 
the usual R-matrix formalism, the scattering (S) matrix is factorized into a hard-sphere 
(hs) factor and an R-matrix factor, the behavior of the latter around resonances defining 
physical parameters: 

U(E) = e 2 ^> = UUE) x \Z%wlB]R{Ey L W = + tP{E), (2) 

where 5 is the phase shift, and S and P are the shift and penetration factors, respec- 
tively Pp. Here, we want to define physical parameters with respect to the potential S 
matrix rather than to the hard-sphere one. By inserting Eq. ^ into Eq. El one gets a 
different factorization for the S-matrix: 

U(E) = U (E) x 1 ~ {[L * (E) - ' R ^ E ^~ l (3) 
y> PoU ^ l-liLW-Bj-i-RvotiE)}- 1 ^^)' 

Comparison of the last factors of Eqs. |2] and El shows that existing R-matrix methods 
transforming physical parameters into formal ones can directly be applied to the hybrid 
case, provided [L(E) — B] is replaced by {[L(E) — B}' 1 — R pot (E)} . In refs. |4|5j . such 
methods are proposed that take into account the interference of the different R-matrix 
poles with each other. Although these methods could be generalized to the hybrid case, 
the interference effects are much less important here since no wide pole is used to describe 
the background. Up to now, I have found sufficient to use the one-pole approximation for 
which, for instance, the physical to formal energy conversion formula reads 

E _ E 2 S(E X ) -B- R pot (E x )\L(E x ) - ij| 2 

A " " l-2R pot (E x )[S(E x )-B]+Rl ot (E x )\L(E x )-By {) 
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Figure 1. Fit of experimental C + a S- 
wave elastic-scattering phase shifts |b|7j by a 
hybrid model consisting of a soft-sphere po- 
tential and one R-matrix pole. 



Figure 2. Soft-sphere C + a effective 
potential (£ = 0) used for the fit of Fig. [TJ 
compared with the pure Coulomb poten- 
tial. 



The conversion formula for the width is too long to be given here. These formulas are 
used for the 12 C + a one-pole example of Fig.^ for which I have checked that the best-fit 
physical parameters of the resonance (Ei = 4.89271 MeV, r\ = 2P(E\)^\ = 1.32 keV) 
are practically independent of the R-matrix radius (above 12 fm, see next section) and 
hence very convenient fitting parameters, contrary to formal parameters. 

3. POTENTIAL CONSTRUCTION BY INVERSION 

An inversion technique based on supersymmetric quantum mechanics has been pro- 
posed in ref. [§] and has been successfully applied to the 12 C + a system in ref. jHj. For 
a given partial wave, this method allows the deductive construction of a potential repro- 
ducing elastic-scattering phase shifts. This method is based on a rational expansion in 
the complex wave-number /c-plane of the partial-wave S matrix 

M K 4- k 

U pot (E = k 2 ) = U (k) J] (5) 

m=l K ™ fc 

where the poles f*t>?ji £1X6 fitted to the phase shifts. They can be determined either by 
resolution of the linear system of equations giving the coefficients of the polynomials 
appearing in the S matrix [9 or with a general minimization procedure ^U]. This last 
method, though less deductive, allows one to easily constrain the poles to stay on the 
imaginary axis, which avoids oscillations in the obtained potential. This effective potential 
is constructed iteratively by algebraic (supersymmetric) transformations of the radial 
Schrodinger equation Vo(r) — > Vj.(r) — > . . . — > Vm{t), with a number of transformations 
equal to the number M of poles of the S matrix. 

The potential Vq and its corresponding S matrix Uq are references which can be chosen 
arbitrarily. For a system without Coulomb interaction, one generally chooses Vq = 0; 
the supersymmetric transformations can then be performed analytically and the final 
potential has a compact expression in terms of exponential functions. For a system with 
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Coulomb interaction, one generally chooses Vq as the pure Coulomb potential. However, 
this has the drawback that individual transformations then create long-range r~ 2 terms 
in the potential. These can be approximately canceled by successive transformations, 
provided most S-matrix poles lie roughly on a circle around the origin of the complex plane, 
but this requires a rather large number of poles and generally creates small oscillations 
in the potential [9,. For instance, to perfectly reproduce the 12 C + a S- wave phase 
shits represented in Fig. [H 10 poles are needed for a pure Coulomb V . This makes 
the construction of the potential rather heavy from the numerical point of view, which 
is not very suitable for applications; moreover, the potential presents slight oscillations. 
An alternative method used here is to construct by supersymmetric transformations the 
nuclear part of the potential instead of the total effective potential. One starts from 
V Q = 0, constructs the nuclear potential Vm, adds Coulomb and centrifugal terms to it 
and only then calculates (nuclear) phase shifts. The fit can only be made by automatic 
minimization procedure in this case. This method presents several advantages: the nuclear 
potential is short-ranged, it has a compact analytical expression and the number of poles 
is small, which facilitates numerical fits. For instance, the fit presented in Fig. ^ has been 
obtained with only two poles at wave number Ki = k 2 = 0.4463? fm -1 . No significantly 
better fit could be obtained without introducing small oscillations in the potential. 

The corresponding effective potential is represented in Fig. |21 together with the pure 
Coulomb potential. V 2 is a purely repulsive soft-sphere potential (no bound or resonant 
state) becoming negligible above 12 fm. This large range is rather surprising: in usual 
R-matrix calculations |7|llj . the range of the nuclear interaction is generally estimated 
between 5.5 and 6.5 fm. Here, a larger R-matrix radius has to be used. A range of 12 
fm could indicate a polarization phenomenon. More probably, it could just mean that 
the soft-sphere potential, which generalizes the hard-sphere potential of R-matrix fits, 
has no strong physical meaning either. This seems confirmed by the fact that I could 
not find, up to now, an ^-independent soft-sphere potential describing the background in 
several partial waves simultaneously. A more promising way is to use a deep potential 
with bound and resonant states corresponding to physical and Pauli-forbidden 16 O states. 
Such potentials can also be constructed by inversion jH] and I plan to test their interest 
for the hybrid model in the near future. 
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